Efficiency of isothermal molecular machines at maximum power 
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FIG. 1. Schematic free energy potential Uo(x) for a two-state 
molecular engine described by a reaction coordinate x under 
the net load force F = Fi - F 2 > 0. 



According to thermodynamics, different forms of work 
can be transformed into one another, with an efficiency 
of at most 100% fl. This lossless limit is achieved with 
a reversible process, i.e., an infinitely slow process. The 
corresponding power output is therefore zero and thus of 
limited interest from a practical standpoint. One of the 
early discussions about efficiency at finite power is at- 
tributed to Moritz von Jacobi around 1840. He realized 
that the output power of an electrical device operating in 
the linear response regime is maximum when the internal 
and external resistors are the same, yielding an efficiency 
of 50%. The Jacobi theorem can easily be reproduced 
in the much more general context of linear irreversible 
thermodynamics: in any engine operating in the linear 
response regime, maximum power is achieved when the 
loading force is equal to half of the stopping force; the 
corresponding efficiency (output power over input power) 
is equal to 1/2. A similar result has been proven for 
the transformation of heat into work, where the maxi- 
mum efficiency, the Carnot efficiency, is again achieved 
under reversible operation, with zero power output. In 
the regime of linear response the efficiency at maximum 
power is again 50% of the Carnot efficiency [2|. More 
recently, in this latter case various explicit results, in- 
cluding bounds for efficiency at maximum power, have 
also been obtained in the nonlinear regime [3l-lll| . In the 



present paper we show that similar results can be derived 
for isothermal molecular machines 12 , . 



Generic model for a molecular motor. We first 
consider a generic model for a molecular motor, namely 
a two-state machine operating along a one-dimensional 
reaction coordinate, see Fig. [TJ The states correspond 
to two minima of an appropriate free energy landscape. 
While a physical energy landscape is expected to be very 
complicated and high-dimensional, the thermally acti- 
vated transitions between the two states will typically 
follow a preferred pathway that connects these states 
via the lowest lying saddle point, the so-called activated 
state. One can project the motion on this pathway and 
introduce a one-dimensional reaction coordinate x with 
corresponding effective free energy potential Uo(x). The 
two "rest" states of the machine, that is, the minima 
in the absence of external forces, correspond to, say, lo- 
cations x = and x = L. The activated state lies at 
an intermediate position x a — XL, < A < 1. In the 
unperturbed phase there are no net transitions, and the 
states 1 and 2 have the same baseline potential value, 
Uo(0) = Uq(L) = 0. The potential has a maximum 
U a = Uo(XL) at the activated state, whose value is typ- 
ically much larger than the thermal energy /3 _1 = fc^T 
(T being the temperature and ks the Boltzmann con- 
stant). In this rest state, the rates, from 1 to 2 and 
k^ from 2 to 1, are equal and given by an Arrhenius law, 
kQ = kg = ko = k exp(— fiUa) ■ We assume a constant 
pre-exponential factor k. 

In the operational regime, that is, in the presence of 
external forces, states 1 and 2 can be identified as "fuel" 
(or "reactant") and "product" states, respectively. To 
transform fuel into product, the machine is subject to 
a driving force F\ which allows it to overcome an op- 
posing but weaker loading force — F2, F2 < F\. These 
forces can be of various physical origins, including chem- 
ical (differences in chemical potentials), electrical (inter- 
nal or external electric fields) or mechanical (e.g., optical 
tweezers, atomic force microscope or optical rotational 
torque). The combined effect of driving and loading is 
a tilting of the potential towards the product state 2, 
U(x) = U (x) - Fx, with F = F 1 - F 2 > 0. In a tran- 
sition from state 1 to state 2, a (scaled) input energy 
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€\ = (3FiL is transformed into a (scaled) output energy 
£2 = fiFiL. The efficiency of this transformation is given 
by 
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Its maximum value, r\ = 1, is reached when the loading 
force i<2 approaches the driving force F\ , and the transi- 
tion from 1 to 2 becomes infinitely slow. In this reversible 
lossless limit the power vanishes. 

In the case of finite and in particular of maximum 
power, the location of the activated state plays a crucial 
role. For a so-called Eyring-like process [12j the activated 
state is very close to the fuel state 1, i.e., A is close to 
zero. The perturbation —Fx barely affects the height 
of the activation barrier that needs to be crossed to go 
from state 1 to 2. The rate also remains essentially un- 
affected, k + ~ fco . However, a maximum barrier increase 
of FL occurs for the backward transition, resulting in 
a rate k~ ps ko exp(— (3FL) (assuming FL <C U a ). On 
the other hand, in the Kramers- like scenario A ~ 1 
k + koexp(f3FL), while fc~ « ko remains essentially 
unaffected. More generally, for a barrier at x a — XL, one 
has k + = fcoexp(Ae) and fc~ = /«oexp[— (1 — -M e L where 
e = ex — £2 = j3FL is the net energy loss or "net load." 
This is proper thermodynamic force (net force divided by 
the temperature) that appears in the entropy production 
and is thus a measure of the distance from equilibrium 

With these explicit expressions for the rates, we turn 
to the output power II given by the output energy €2 
multiplied by its net rate of production, II = fce2, with 
k = k+-k- = k(e) = k [e Ae - e^ 1 "^] . To specify the 
condition of maximum power we set dll/de2 = 0, which 
yields the unique solution 
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A(l - e- £ ) + e- e 
= e +(i-A)e 2 + (i-A + AV + 0( e 4 ). (3) 

This result in Eq. (JlJ yields one of the central results of 
this paper, namely the efficiency at maximum power: 
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+ 0(£ 3 ). (5) 



We point to a number of revealing observations. The 
first term of the expansion |5]) is the prediction of linear 
irreversible thermodynamics of efficiency at maximum 
power equal to 1/2. The associated relation between the 
forces, 2^2 = Fi is obtained from the first term in ex- 
pansion ([3]), e 2 = e = ei - £2- 

Turning to the next order corrections in Eqs. Q and 
(JSJ), the coefficients vanish in the symmetric case A = 1/2, 



reminiscent of a similar property for thermal machines 
fTol j . Note also that the coefficient of the term propor- 
tional to e in (JSJ goes from a maximum value 1/8 at 
A = to the minimum value —1/8 for A = 1, switching 
from positive to negative values at A = 1/2, again remi- 
niscent of an analogous feature in thermal machines 11 1 . 
The first two terms of the expansion ([5J were also derived 
in [l4j |. but the connection with the physically relevant 
parameter A was not made (see, however, [Hj]). 

The efficiency 77* at maximum power is a function of 
A and e. One easily verifies that 77* is a monotonically 
decreasing function of A for given e > 0. The upper 
limit is the efficiency 77^ = 77* (0, e) < 1 of the extreme 
Eyring-like scenario, and the lower limit is the efficiency 
< 77* (1, e) = rf of the extreme Kramers- like case: 
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= Ve- (6) 



We next consider the e dependence, starting with the 
variation of the bounds. The Eyring-like efficiency 77}^ 
increases monotonically from 1/2 when e — > to rj^ = 1 
when e — > 00. The Kramers- like efficiency 77 * K decreases 
monotonically from 1/2 when e — ¥ to r]* K = when 
e — > 00. The variation of 77* between these bounds de- 
pends on A. When A > 1/2, 77* decreases monotonically 
from 1/2 when e = to when e — > 00. The system is in 
the product regime (Kramers- like), and 77* behaves much 
like the Kramers- like limit, never rising above the lin- 
ear response value 1/2 [dotted curve, short-dashed curve, 
and filled circles in Fig. [2ja)]. On the other hand, when 

< A < 1/2 the system is in the fuel regime (Eyring-like) , 
and 77* starts at 1/2 when e = 0, rises to a unique max- 
imum, and then decreases to as e —> 00 [long-dashed 
curve, solid curve, and open circles in Fig.[^a)]. The "op- 
timal" value of the efficiency at maximum power occurs 
at the net load value e which solves the transcendental 
equation obtained by setting the derivative of (j4]) with 
respect to e equal to zero, (1 — e e ) [A(l — e c ) — 1] = ee e . 
Each point along the curve in the inset of Fig. [D^a) is 
associated with a different value of e. High efficiencies at 
maximum power require the system to operate very near 
the fuel state. Thus, for instance, referring to the figure, 
the maximum of the A = 0.1 curve (solid) is 77* = 0.69 
and occurs when the net load is e = 3.19. A maximum 
efficiency of say 77* = 0.9 requires that the net load be 
e = 5.68 and that the motor operate at A = 0.016. 

We can repeat our analysis for e < 0, with net transi- 
tions going from state 2 to state 1. Indeed, many motors, 
including ATPase, can operate in reverse. The inter- 
change of 1 and 2 corresponds to a replacement of A by 

1 — A. Hence the above theory indicates that, when con- 
sidering both modes of operation, at least one of them 
has A > 1/2, with a corresponding efficiency at maxi- 
mum power less than 50%. For A = 1/2 the engine works 
equally well at maximum power in forward and reverse 
modes, at 50% efficiency. 
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FIG. 2. Efficiency at maximum power, rf , as a function 
of net load e. In all panels the open circles represent if E 
and the filled circles rf K . Panel (a): two state model [cf. 
Eq. QJ] with A = (long dashes), 0.1 (solid), 0.5 (short- 
dashes), and 0.9 (dots). The solid and long-dashed curves 
exhibit a maximum, clearly seen on this scale in the solid 
curve, while the two lower curves are monotonic. The inset 
shows the maximum value of rf as a function of A, with < 
A < 1/2. The two panels in (b) are for the periodic square 
well potential, and those in (c) are for the sawtooth potential. 
For the former 5 — 0.1 in the left panel, 0.5 in the right panel. 
In both panels, rf is shown for three barrier heights: f3U a = 1 
(solid), 3 (dashes), and 9 (dots). For the sawtooth potential 
the curves in the left panel are all for fiU a = 10 and in the 
right panel for (3U a = 100. The four curves are for different 
values of the potential asymmetry parameter A: A = (long 
dashes), 0.1 (solid), 0.5 (short dashes), and 0.9 (dots). 



Generalized model. The above generic model as- 
sumes exponentially difficult crossing of the transition 
state. With a more general analysis, we suppose that the 
motion projected on reaction coordinate x can be de- 
scribed as an overdamped one-dimensional diffusion pro- 
cess in a potential U(x) = Uq(x) — Fx: 
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dt 



dx 



(7) 



Here 7 is the viscous friction coefficient and £ is Gaus- 
sian white noise, (£(t)) = 0, and (£(*)£(*')) = - 
The diffusion coefficient is related to the viscous friction 
coefficient and the temperature by the Einstein relation 



D = ksT/^. This description relates transition rates 
to more fundamental parameters than the earlier phe- 
nomenology. It has been shown to be a very useful tool to 
describe the response of molecular machines [l6| , and it 
has the further advantage of being analytically tractable. 
It also does not require extremely high activation barri- 
ers. Indeed, when the potential barriers are comparable 
to the thermal energy fee T, it is no longer appropriate to 
identify the minima of the potential (formerly x = and 
x = L) as states between which transitions take place. 
We therefore replace the two-state scenario by an ex- 
tended coordinate x £ [—00,00], and consider a baseline 
potential Uq(x) periodic in x with period L. 

Upon application of a driving force F% and a load force 
F2, the net force F = F\ — F2 > induces a steady 
state current in the tilted potential U{x) = Uq{x) — Fx 
with average velocity V — (dx / dt) along the positive in- 
coordinate. The transformation of driving energy e\ — 
pF\L per period into loading energy €2 — ftFiL per pe- 
riod takes place at a net rate k = V/L, with output power 
IT = ke2 = Ve2/L. The efficiency of the transformation is 
again given by 62 jt\ , see Eq. ([T]) . We note in passing that, 
besides being a natural model for cyclic molecular motors 
such as ATPase, overdamped Brownian motion in a tilted 
periodic potential also provides a relevant description in a 
large number of other physical situations [17], including 
Josephson junctions, rotating dipolcs in external fields, 
particle separation by electrophoresis, transport in tubes 
of varying cross-section, and biophysical processes such 
as neural activity and intracellular transport. 

Turning to the issue of efficiency at maximum power, 
we first derive results that hold for arbitrary potential. 
We suppose that the average velocity V can be written as 
a power series in F. Since V vanishes for F = 0, there is 
no constant term in the expansion. For comparison with 
the previous results, we consider the rate k = V/L for 
moving over one period L, and write the corresponding 
power series in terms of e = f3FL, k = k(e) = V/L = 
die + <22£ 2 + a^e 3 + 0(e 4 ). Maximization of the output 
power IT = fc(e)e2 with respect to the loading energy 
£2 gives the following expansion for the output yield at 
maximum power [compare with Eq. ([3])]: 
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The corresponding efficiency reads [see Eq. ||5J]: 
1 a? (3a2 — 4aia3) 
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This expansion features the familiar 50% efficiency in 
the regime of linear response. Turning to the nonlin- 
ear regime, we note that just as in the generic two state 
model, the next order correction vanishes (02 = 0) when- 
ever the system has left-right symmetry for the velocity, 
V(F) = -V(-F). 
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The average steady state velocity for overdamped mo- 
tion in a tilted periodic potential is given by |l8[ 



v(t) . L L MP-O . 



(10) 



This expression in principle allows us to find the power 
II = Ve2/L, and hence makes it possible to explicitly 
identify the regime of maximum power and its corre- 
sponding efficiency. 

Firstly, we rederive the results for the two-state model 
with a potential Uq with a dominant high maximum in 
each period, say at x = x a = XL (modulo L), and a 
unique minimum at x = (modulo L). The dominant 
contribution to the double integral in Eq. (fT0| comes from 
the region around the (x,y) point for which Uq(x + y) 
reaches a maximum and Uq(x) a minimum. This point 
lies at x + y = x a = XL and x — 0, and consequently y = 
XL. The e-dependence of the denominator is therefore 
of the form exp(— Ae). It then follows directly that V ~ 
k ~ {1 — cxp(— e)} exp(Ae), and the resulting power is 
identical to that for the two-state model. 

Secondly, we note that universal conclusions can be 
drawn using Eq. (|10|) even without an explicit evalua- 
tion of the integrals. In particular, we can identify the 
coefficients aj in Eq. ([3]), 
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where 
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Since < 1% < Jo, the coefficient a 2 /4ai = (Jl/Jo — 
l/2)/4 of the linear term lies between —1/8 and 1/8, as 
was the case for the 2-state model, see Eq. ([3]). Further- 
more, the coefficient is zero for a potential with left-right 
symmetry, that is, when there exists a point xq for which 
U (x - x ) = Uq(x q - x). 

To proceed further, it is in general necessary to in- 
voke numerical calculations because the integrals in (|10p 
can not be performed analytically for a general potential 
Uq{x). However, analytic results can be obtained in some 
limits (such as the dominant high maximum case consid- 
ered above), or for some specific shapes of the poten- 
tial. Examples include the square well, Uq(x) = U a , x 6 
[0,<5J] modulo L, and Uq(x) = 0, x <E [5L,L] modulo J, 
and a saw-tooth potential. To illustrate, we quote the ef- 



ficiency at maximum power for the square well potential, 
- l) 

(13) 
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(1 - e - e £ ) + c [l - Sd 1 -^ + (6 - l)e 5e ] ' 
where c = 4sinh (/3U a /2). The efficiency rf for 
the square-well and sawtooth potentials are shown in 
Figs, (lb) and (c). 

Finally, we mention a generic behavior for any poten- 
tial with finite maxima. For very large driving, the bar- 
rier(s) of the potential Uq, indeed the entire potential, 
become irrelevant and one goes back to a linear model 
with V — F/"f = ait with a\ = D/L, implying that 
rf = 1/2. This return to the linear scenario can be seen 
in some cases in Fig. [2j in others one must go to higher 
values of e than those shown in the figure. 

Closing perspective. Is the issue discussed in this 
paper, maximizing power with respect to the load, a rele- 
vant criterion in practice? In the case of thermal motors, 
this seems to play a role, at least from an engineering 
point of view, since power plants operate under condi- 
tions in general agreement with this criterion (llj . We 
hope that the present paper will lead to a re-examination, 
from the perspective of maximum power, of the much 
larger class of isothermal engines, including the impor- 
tant class of molecular motors. 
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